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Abstract 

We study the thermodynamic limit of random partition models 
for the instanton sum of 4D and 5D supersymmetric U(l) gauge theo- 
ries deformed by some physical observables. The physical observables 
correspond to external potentials in the statistical model. The parti- 
tion function is reformulated in terms of the density function of Maya 
diagrams. The thermodynamic limit is governed by a limit shape of 
Young diagrams associated with dominant terms in the partition func- 
tion. The limit shape is characterized by a variational problem, which 
is further converted to a scalar- valued Riemann-Hilbert problem. This 
Riemann-Hilbert problem is solved with the aid of a complex curve, 
which may be thought of as the Seiberg-Witten curve of the deformed 
U(l) gauge theory. This solution of the Riemann-Hilbert problem is 
identified with a special solution of the dispersionless Toda hierarchy 
that satisfies a pair of generalized string equations. The generalized 
string equations for the 5D gauge theory are shown to be related to 
hidden symmetries of the statistical model. The prepotential and the 
Seiberg-Witten differential are also considered. 
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1 Introduction 



This is a sequel to our previous paper pQ on the correlation functions of loop 
operators of the 5D supersymmetric 17(1) gauge theory in the Q background 
on R 4 x S 1 [2J. As shown therein, a generating function of those correla- 
tion functions coincides with the partition function of the so called "melting 
crystal model" deformed by a set of external potentials. The melting crystal 
model is a statistical model of random plane partitions (3D Young diagrams). 
Actually, by the technique of diagonal slicing [3], the partition function can 
be converted to a sum over ordinary partitions (Young diagrams) . The latter 
may be thought of as a g-analogue of the instanton sum (Nekrasov function) 
of 4D M = 2 supersymmetric gauge theories [I]. The undeformed melting 
crystal model was first proposed as a statistical model of the amplitude of 
A-model topological strings [5J. A slightly modified model was studied as 
the instanton sum of 5D supersymmetric gauge theories [HI [7j. The fully 
deformed model was introduced from the point of view of integrable struc- 
tures [S] (see also the review [H]). It was, however, not clear what physical 
meaning the external potentials have in the context of the 5D gauge theory. 
Our previous paper [1] presented an answer to this question. 

Another issue addressed therein is the thermodynamic limit of the melt- 
ing crystal model, from which we attempted to derive a 5D analogue of the 
Seiberg-Witten theory [10J. In the case of 4D M = 2 supersymmetric gauge 
theories, such a "microscopic" derivation was first achieved by Nekrasov and 
Okounkov [11] . and later extended to deformed models by Nekrasov and Mar- 
shakov [121 IT3] . A clue to this problem is the notion of "limit shape" . In the 
thermodynamic limit, the partition function is dominated by Young diagrams 
of macroscopic sizes and a particular (rescaled) shape. This limit shape can 
be determined by a variational problem. The variational problem, in turn, 
can be converted to a Riemann-Hilbert problem for a complex analytic func- 
tion. Solving the Riemann-Hilbert problem, Nekrasov and Marshakov [T2|[T3] 
could derive the Seiberg-Witten theory in the presence of external potentials 
(which they identified with local observables Tr <p k+l / {k + 1), k = 1, 2, • • • , 
in the M = 2 chiral multiplet). Inspired by their work, we studied the 
thermodynamic limit of the deformed melting crystal model, using a slightly 
different method that was developed for the undeformed model [3 HI]- Unfor- 
tunately, our calculations for the deformed model were wrong; we considered 
an inappropriate Riemann-Hilbert problem, which led to wrong results. 

Our primary problem in this paper is, therefore, to formulate a correct 
Riemann-Hilbert problem for the deformed melting crystal model. We ad- 
dress this problem in a wider context. Namely, we treat both the 4D and 
5D theories by the same method, thereby hoping to compare these two cases 
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on a common ground. Moreover, this enables us, simultaneously, to com- 
pare our method with the method of Nekrasov and Marshakov. As it turns 
out, one can formulate (and solve) an appropriate Riemann-Hilbert prob- 
lem for a function W(z) that amounts to the so called "loop amplitude" or 
"resolvent function" in the theory of random matrices. The solution is ac- 
companied with a complex curve (a double covering of the z-plane), which 
may be thought of as the Seiberg-Witten curve of the U(l) gauge theory. In 
the case of the deformed 4D theory, this result fully agrees with the result of 
Nekrasov and Marshakov. 

We further present an interpretation of these calculations in the language 
of the dispersionless Toda hierarchy [T5l [T6] . The dispersionless Toda hier- 
archy is a long-wave limit of the Toda hierarchy [T7] (see the review [TS] for 
more details). It is known [TT | IT3" j l8] that both the 4D and 5D Nekrasov func- 
tions are related to the ID Toda hierarchy (higher time evolutions of the Toda 
lattice), hence to the Toda hierarchy (higher time evolutions of the 2D Toda 
fields) as well. On the other hand, as already observed by Marshakov and 
Nekrasov (loc. cit.), the thermodynamic limit corresponds to the long- wave 
limit of the Toda lattice, and the solution of the aforementioned Riemann- 
Hilbert problem in the 4D case gives a solution of the dispersionless ID Toda 
hierarchy. We consider this issue from a different point of view, namely, the 
notion of "generalized string equations" in the Toda and dispersionless Toda 
hierarchies. This notion stems from string theories [T3j I2TH |2~T| |2"2~| [25] , and 
found applications in interface dynamics [2U ESI EE] , which eventually led to 
a detailed description of a quite general class of solutions of the dispersion- 
less Toda hierarchy [27] • Employing this notion and techniques developed in 
earlier studies, we show that the solutions of the Riemann-Hilbert problem 
for both 4D and 5D cases can be (almost uniquely) characterized by a pair 
of generalized string equations. "Almost" means that the generalized string 
equations for the 5D case are suffered from a kind of ambiguity. The meaning 
of this ambiguity is left for a future study. 

Let us mention that the notion of generalized string equations is also use- 
ful for a statistical model of the double Hurwitz numbers [28]. Remarkably, 
this statistical model and the melting crystal model have a common alge- 
braic structure in their hidden symmetries. Following this analogy, one an 
derive the generalized string equations of the 5D theory from those hidden 
symmetries of the partition functions. We discuss this issue in the end of the 
paper. 

This paper is organized as follows. Sections 2, 3 and 4 show the formu- 
lation of the random partition models, the prescription of thermodynamic 
limit and the derivation of Riemann-Hilbert problems. In Section 2, the de- 
formed random partition models are formulated. A 2D free fermion system 
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is used to formulate the main statistical weights and the external potentials 
in a compact form. In Section 3, the partition functions are reformulated in 
terms of the density function of Maya diagrams. The energy functionals for 
the density function are specified. In Section 4, the thermodynamic limit is 
formulated and the asymptotic form of the energy functionals is determined. 
The variational equation for the limit shape is derived and converted to a 
Riemann-Hilbert problem. Sections 5, 6, 7 and 8 present the main results 
on solutions of the Riemann-Hilbert problem and their relation to the dis- 
persionless Toda hierarchy. The solution for the 4D theory is constructed in 
Section 5, and examined in Section 6 from the point of view of generalized 
string equations. The solution for the 5D theory is considered in the same 
way in Sections 7 and 8. Section 8 also presents a derivation of the gen- 
eralized string equations of from hidden symmetries of the melting crystal 
model. Section 9 is devoted to concluding remarks on the prepotential and 
the Seiberg-Witten differential. 



2 Random partition models of U (1) gauge the- 
ories 

2.1 Fermions 

We use the same formulation of fermions as our previous papers [H [9j Q] . Let 
ip n ,ip*, w 6 Z, denote the Fourier modes of 2D complex fermion fields 

neZ neZ 

They satisfy anti-commutation relations 

The charge-s sector of the Fock space has the ground states 

(s\ = (-Oo| ■ ■ • C-lVC: I s ) = Ip-slp-s+l • • • | - OO) 

and excited states 



(ji, S\ = (-00| • ■ ■ V£ 2+s -lV£ 1+s , \H, S) = Ip-^s^-^-s+l • • • I - 00} 

labelled by the set V of all partitions fi = fi\ > Hi > ■ ■ ■ > 0, of 

arbitrary lengths. 
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Let us introduce the fermion bilinears 



neZ 



neZ neZ 



Jfc's span a U(l) current algebra, and L and W are zero-modes of Vira- 
soro and algebras. J , L and Wo are diagonal with respect to the 

orthonormal basis (//, s|, \(J>,s). The diagonal elements can be calculated as 



(/!, S\j \fi, S) = S, 

s ■ 

( 

s(s + l)(2s + 1) 



(/i,s|L |/i,s) = | A*l + g ' (2- 1 ) 



(//,s|W |//,s) = k(a*) + (2p+l)|//| + 
where 

Ia*i = 5^^' k ^ = ~ 2i + 1 )- 

z>l i>l 

We use this fermionic language to formulate the random partition models 
of instanton sums for the 4D and 5D U(l) gauge theories. The models are 
deformed by an infinite number of external potentials $/c(/i, s), k — 1, 2, • • • . 
The partition function of each model there becomes a function Z s (t) of the 
coupling constants t = (ti,t 2 , ■ ■ ■ ). The coupling constants play the role of 
time variables in an underlying integrable hierarchy. 

2.2 Partition function for 4D theory 

We start from the the fermionic formula 

Zf>(t) = (s\e R - 1J iAl Lo e H We h - 1J -i\s), (2.2) 
where A and h are positive constants, and H(t) is a linear combination 



H{i) = Y J hH k 



k=i 

of the zero-modes 

h k+i 



Hi. = 



T j2( nk+1 -( n -^ k+1 y^n- 



k + _ 

neZ 
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of (a slightly modified version of) the ordinary algebra. 

Since the action of e Jl , K^ L °e H ^ and e J_1 preserves the charge- s sector, 
one can insert the partition of unity 

1 = \fjL, s\ (in charge-s sector) 

between these operators, where V denotes the set of all partitions of arbitrary 
length. This leads to an expansion of Zf^{t) into a sum over partitions. 
Moreover, since A 2L °e H w is diagonal, this becomes a single (rather than 
double) sum: 

Zfit) = ^(3|e^ 1Jl |/i, S )(/i, S |A^e^V^)(/i^|e /8 " 1J ' 1 k). 
i-tev 

Let us specify the building blocks of this sum. Firstly, e J±1 are known to 
act on the ground states as 

where dim /x denotes the degree dim of the irreducible representation 
of the symmetric group S r , r = \fi\ [29]. dim/i has the hook length formula 

dim^= YlKD)' 1 , (2.4) 

□ 6/1 

where h(Di) denotes the hook length of the cell □ in the Young diagram of 
/j. Secondly, the diagonal elements of A 2L °e H ^ can be factorized as 

(fi,s\Al Lo e H V\n,s) = a%I+^+V^, (2.5) 

where 

oo 

$(*,//, s) = yjk^kdi^), $ fe (/i,s) = ((J,, s\H k \n, s). 

k=l 

$fc(/i, s) has the formal (apparently divergent) expression 
fcfc+i 00 

s) = j-^j J2 H s + tM - i + - (« + V* ~ *) k+1 ) > (2-6) 
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which can be reorganized to a finite sum as 
fcfc+i 00 

i=l 

i=l 

Thus Z^ D (t) can be converted to the sum 

Zf{t) = 2 Ajf l+s(s+ (2.8) 

over all partitions. This is exactly the Nekrasov function of the 4D M = 2 
supersymmetric U(l) gauge theory deformed by external potentials 



2.3 Partition function for 5D theory 

We define a 5D analogue of ( 12.21) as 

Z™(t) = (s\G + Q Lo e H ^G_\s), (2.9) 

where q and Q are constant with < q < 1 and < Q < 1, is a linear 
combination 

00 

H(t) = J2tkH k 

k=i 

of the zero-modes 



H k = J2q kn :ij_ n r n . 



neZ 

of the quantum-torus algebra [9], and G± are the transfer operators 

/ 00 fe/2 



J^*(l-g*) 

of Okounkov and Reshetikhin [3]. 

Inserting the partition of unity between operators, one can expand Z^ D (t) 
into a sum over all partitions. The building blocks of this expansion can be 
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calculated as follows. According to Okounkov and Reshetikhin [3], G± act 
on the ground states (s\, \s) as 

(s\G + = Y t (»,s\s (t (q»), G.\s) = J2^(Q P )\^s), (2.10) 

where s M (g p ) is the special value of the Schur function s^Xi, X2, ■ • ■ ) of infinite 
variables [30] at 

This special value has the hook length formula 

S^qP) = J](g-MD)/2 _ /W)-^ ( 2 .11) 

Den 

which may be thought of as a g-deformation of the hook length formula (12. 4p 
of dim jj,. The operator Q L °e H ^ in the middle is diagonal, and the diagonal 
elements can be factorized as 

(ji, s\Q Lo e Hit) \^ s) = QM+^+m e m^)^ (212) 

where 

oo 

$(*,//, s) = y]tk®k(l*,s), ®k(n,s) = (fi,s\H k \n,s). 

k=l 

s) has the formal expression (apparently valid for \q\ > 1) 

oo oo 

$ fc (^ s) = J2 q k(s+ ^- i+1) - Q ki ~ i+1) , (2-13) 
which can be reorganized to a finite sum as 

OO ^/^ jTj^x 

Sfcfo s) = J2(q k{s+ ^ l+1) - + q ; _ ! J - (2.14) 

i=i ^ 

Thus Zg D (t) can be expanded to a sum over all partitions as 

Zf{t) = s^q p ) 2 Q^ +s{s+1)/2 e^ s) . (2.15) 

This is a 5D analogue of the deformed 4D U(l) Nekrasov function (12. 8p . 
Although one can further insert a Chern-Simons term [6], we consider this 
simplest model in this paper. 
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3 Energy functional for density function 

3.1 Density function of Maya diagram 

We now recall the notion of density function of Maya diagrams [7J [H] that 
play a central role in our approach to the thermodynamic limit. 

The Maya diagram of \i is a configuration of particles located at the points 
x = Hi — i, i = 1,2, on a line. Let p M (x) denote the density function 



PA X ) 



-fii + i) 



8=1 



of this particle configuration. The backward difference 

&Pn(x) = p^x) - p^x - 1) 

becomes, up to a constant factor, the second derivative of the so called "pro- 
file function" f^x) [TT] of the (45 degrees rotated) Young diagram: 

Ap,(x) = -\f;{x). (3.1) 
The first few moments of Ap^(x) can be explicitly calculated as 

oo 

dxAp^x) = -1, 

■oo 
oo 



dxxAp^(x) = 0, 
dxx 2 Ap fl (s) = —2\n\, 

oo 
oo 

dxx 3 Ap^(s) = —3k(p). 



(3.2) 



The density function p^ s (x) for the charged partition (fi,s) is a parallel 
transform of Pfj,(x): 



■oo 



PnA x ) = Pm( x - s ) = y~]8(x- s - iij + i). 



i=l 
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The first few moments read 



/oo 
dxAp^ s (x) = -1, 
-oo 

/oo 
dxxAp^ s (x) = -s, 
-oo 

/oo 
dxx 2 Ap^ s (s) = -2|//| - s 2 , 
-oo 

/oo 
dxx 3 Ap^ s (s) = —3k(p) — 6s\p,\ — s 3 
-oo 



(3.3) 



3.2 External potentials in terms of density function 

Let C^ D (/i, s) denote the k-th moment of — Ap^ s (x): 

/oo 
dxx k Ap^(x). 
-oo 

These moments have a formal expression 

oo oo 

0*V *) = j> + ^ - i + 1)" - j> + ^ - 0*. (3-4) 

i=l i=l 

which can be converted to a finite sum as 

oo 

0^(^=0* + *-* + !)*-(* -* + !)*) 

i=l 

oo 

-X)(( a + ^-0*- (« -<)*)+«*■ (3-5) 

i=l 

The external potentials 

f-k+l 00 

**V s) = ^ £ ((a + m - i + - ( a + ^ - i) k+1 ) 

i=l 

of the 4D gauge theory can be thereby expressed as 

fU.) = = _£ d ^ Ap „ (l) . ,3.6) 

One can find a similar expression of the external potentials for the 5D 
gauge theory. Consider the g-analogue 



/oo 
-oo 
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of Of D (p, s). The formal expression 

oo oo 
Ol D {fI, S) = J2 - ^ g *(-+W-<) (3.7) 

i=l i=l 

can be converted to a finite sum as 

oo 
i=l 

oo 

- XI (g fc(s+Ml "° - + q ks . (3.8) 

i=l 

The external potentials 



^f(A*,a) = E9* ( * +W " <+l) -E9* 



i7 fc(-i+i) 

i=l i=l 

of the 5D gauge theory are related to these g-moments as 



A; 



fe /-oo fe 

^ Q J — OO * 



(3.9) 



3.3 Hook product in terms of density function 

By the hook length formulae (12.41) and (12. lip , the main part of the Boltzmann 
weights in (12. 8p and fl 2 . 1 5 [) can be rewritten as 

(J^)>^)W>, (3,0, 
s,(q p ) 2 Q M = q~ KM/2 J] <5^ (D) (1 - g MD) )~ 2 - (3.11) 

Let us recall the general formula [lj 

X)/(/i(D)) = ~ IT dxdyg{\x-y\)Ap,(x)Ap,{y), (3.12) 

where /(x) is an arbitrary function and g(x) (referred to as the kernel func- 
tion) is a function that satisfies the conditions 

g(x + l)+g(x-l)-2g(x) = f(x), g(0) = 0. (3.13) 

By this formula, the logarithm of (I3.10P and (13.111) can be converted to 
quadratic functionals of Ap^ s as follows. 
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4D theory [TT]: One first obtains the expression 
dim fi 



where 



. . x(x - 1) . A . „ , . N 
94d(x) = i og _£ + l og G 2 (a; + 1). 



G^O^) is Barnes' G-function, the second member of the multiple Gamma 
functions G n (x),n = 0, 1, • • • . Note that (?4d(^) satisfies the equation 

fix 

g{x + 1) + g(x - 1) - 2g(x) = log — . (3.14) 

A o 

Substituting x — > x — s, one can rewrite Ap M to Ap^ s as 

- log (mj^j A l M = J J dxdyg m (\x - y\)A Pfl , s (x)A Pfl , s (y$.15) 

5D theory [1] : Also using the last formula of (13.21) to the term (k( P )/2) log 
one first obtains the expression 

- \og Sfl (q p ) 2 Q M = dxdyg 5D (\x - y\)Ap^x)A PpL {y) 

J J— oo 

logg 



dxx A Pll {x) 



oc 



6 

where 

x(x - 1) x(x 2 - 1) 
£sd (ar) = log Q — log q 

+ ~ log(l - g) + logG 2 (x + 1; g). 

G2(x;q) is the second member of the multiple g-Gamma functions 
G n (x; g), n = 0, 1, • • • [31]. Note that (7sd(^) is chosen to satisfy the 
equation 

g(x + 1) + g(x - 1) - 2^(x) = - + log(l - (3.16) 

Substituting x — > x — 5, one can rewrite Ap^ to Ap M;S as 

-log s^qPfQ^ = // dxdy# 5D (|:c - y\)A P », s {x)A PfljS (y) 

, °° ,00 (3.17) 
logg p 3 

dx(x - s) J Ap AtiS (x), 
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3.4 Energy functional for density function 

Thus the main part of the Boltzmann weight becomes the exponential of a 
quadratic functional of Ap^ s . The external potentials are already shown to 
have such a functional expression by (13. 6p and (13. 9p . Thus, up to a sim- 
ple factor independent of Ap^ s , the Boltzmann weight can be cast into the 
standard form e~^ PA1 ^ with the following energy functionals S[p\: 

• 4D theory: 

£t D [p] = dxdyg m (\x - y\)Ap(x)Ap(y) 

J J —oo 

/oo / °° 

E ""»-7i 1 Ap(x). (3.18) 
-OO \ 



fc + 1 

v fc=l 



5D theory: 

O] = / / ^^5d(|x - y|)Ap(ar)Ap(y) 

•/ J —oo 

/OO /-CO / OO \ 

dx(x - s) 3 a P (x) -J dx[j2 r^ki kx ) a p( x ) 



log ? 
6 



(3.19) 



We now have the following expression of the partition functions as statistical 
sums over the set V s of all density functions of the form p = p^ s , \i e V: 

Zf>(t)=k s Q {s+1) Y J Z- e " [p \ (3-20) 

pev s 

t k q k 
1-5 , 



Zf{t) = ew[J2f ± - k ) Q s{s+1)/2 £ e "^ W " ( 3 - 21 



4 Thermodynamic limit and Riemann-Hilbert 
problem 

4.1 Thermodynamic limit for 4D theory 

To achieve the thermodynamic limit for the 4D theory, we rescale s and t^s 

as 

s -)> fr^s, tjfe -)> /r 2 t fc , (4.1) 
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and let h — > 0. The partition function is then dominated by partitions p, 
with = 0(tV 2 ). The rescaled profile functions 

frU,K-^s{h~ l u) = hf^tr 1 ^ - s)) 

of these dominant partitions then tend to a continuous "limit shape" /j (u). 

To formulate this limit shape in the language of the density functions of 
Maya diagrams, we assume that the dominant contribution to the partition 
function is due to density functions of the form 

p{tT l u) = p (0) (u) + 0(H) as fi -> 0, (4.2) 

where p(°'(u) is a continuous function, monotonously decreasing, and differ- 
entiable in a weak sense, and the derivative p^'(u) satisfies the constraints 

POO 

dup {0) '(u) = -1, / duup (0) '(u) = -s. (4.3) 



Note that these constraints stem from the moment formula (13.31) . 

Under these assumptions, we can find the leading part (of order fir 2 ) of 
the energy functional S^[p] explicitly. Note that the kernel function g 4D in 
the rescaled coordinate behaves as 

g iD (fi- l u) = fi- 2 g ( £(u) + OifT 1 ), (4.4) 

where 

(o)/ \ u 2 ( u 3\ 

Also note that Ap(x) in the rescaled coordinate can be expressed as 

Ap(fiT l u) = hp {0) '(u) + 0(fi 2 ). (4.5) 

Consequently, the quadratic part of £t°[p] turns out to have the following 
asymptotic form: 

\ \ dxdyg iD (\x - y\)Ap(x)A(y) 

J J —CO 

P POO 

= fi 2 dudvg { ^(\u - v\)p^'(u)p^'(v) + 0(fi- 1 ). 



lr This means that the derivative p^' can be singular at some exceptional points, but 
that the various integrals containing p(°>' are still meaningful. 
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On the other hand, the linear part, after rescaling ty. — > h 2 tk, behaves as 

= ^ 2 du fa ^ j + O^- 1 ). 

Thus the energy functional itself can be expressed as 

^ 4D 2 JP] = ^ 4D( ° ) [P (0) ]+0(^ 1 ), (4.6) 

where 

/» p oo 

C (0) [P (0) ] = / / dudvg^(\u - v\)p^(u)p^'(v) 

J J —oo 



(«)■ (4-7) 



4.2 Riemann-Hilbert problem for 4D theory 

The density function 

pl 0) H = 1 "f )M 

of the limit shape is obtained as a minimizer of £ t 4D(l V )] under the con- 
straints f !4.3p . Following the well known ansatz, we assume that p* (it) is 
constant outside an interval [«o, wi]. This amounts to the "one-cut" condi- 
tion in the large- N limit of random matrix models. Since the test functions 
(it) in the minimization problem are monotonously decreasing continuous 
functions with the obvious boundary conditions 

lim p (0) (w) = 1, lim p (0) (w) = 0, (4.8) 

u— >— 00 u— >+oo 

this means that 

Pf\u)=\ l n (4.9) 
10 tor u > u\. 

If we limit the test functions to those with the last property (14. 9p . the 
variational equation (with Lagrange multiplier v) 

6 Uf m [p (0) ] + v ( [°° duup^\u) + s)) = 
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reduces to solving the integral equation 

Tu {£ dvg{ ^ u - u i^ (0) ») + + \ = {uo - u - (4 - 10) 

where 

oo 
k=l 

under the constraints (14. 3 p and the support condition H 

p(°)'( M ) = for u [uq, U\\. (4.11) 

Let us stress that the endpoints uq, U\ are "dynamical", namely, determined 
along with the minimizer p* (u). 

Actually, rather than solving (I4.10p directly, we consider the equation 

PP f * dvgS'(\u - v\)pP»(v) + ^ = (u < u < « x ) (4.12) 

obtained by differentiating (I4.10p . Note that the integral on the left hand 
side is now interpreted to be the principal part (indicated by the notation 
"PP"), because the kernel function 

(o)///i I \ i \u — v\ 
^4D (\ u - v \) = lo g— r — 

has a logarithmic singularity at v — u. 

We now introduce the complexified kernel function 

(0)/// \ , z-v 
AO 

and construct the complex analytic function 

W{z)= ^ dvgS'(z-v)p^(v) 

from the derivative p* (it) of the minimizer. W(z) is a holomorphic function 
on the cut-plane C \ [u , Ui] with the boundary values 

W(u ± i0) = PP r dvgZ'(\u - v\)p^'(v) T Tripsin) (4.13) 

2 Actually, the derivative p^' (u) of the minimizer turns out to have singularities at the 
endpoints uq,ui of the support. The integrals in (14. 10)) and (|4.3p . however, have finite 
values. 
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along the real axis. One can thereby rewrite (14. as 

W(u + iO) + W(u - iO) = -V'(u) (m <m<mi), (4.14) 
and supplement this equation with the equation 

W (u + zO) - W(u - iO) = | ° ^> Ul ) (415 ) 

I — 2ni (u < u ) 

off the support of p^' . Moreover, from the asymptotic form 

94d"( z ~ v ) = 1°§ 7 — I" 0( z ~ 2 ) as z — )• oo 

A 2 

of the kernel function and the constraints ( 14. 3p . one can see that W(z) should 
satisfy the boundary condition 

W(z) = -log-^- + - + 0{z~ 2 ) as^^oo. (4.16) 
A z 

Conversely, the constraints ( 14.31) can be recovered from this boundary con- 
dition. 

(I4.14p . (14.151) and (I4.16P altogether form a kind of Riemann-Hilbert prob- 
lem. Once a solution of this Riemann-Hilbert problem is found, the minimizer 
pi '* can be obtained from the imaginary part of W(u ± 0) as 

P ?Hu) = T ImW{U±i0) . (4.17) 

7T 

Let us mention that one can translate the variational equation ( I4.10p 
itself, too, to a Riemann-Hilbert problem. This Riemann-Hilbert problem is 
formulated for the primitive function El 

G(z) = J dzW(z) 

(with suitable normalization of the integration constant) of W(z). One can 
solve this Riemann-Hilbert problem in much the same way as that of W(z). 
Actually, this is exactly how Marshakov and Nekrasov treated the thermo- 
dynamic limit [T21 [T5] . Speaking more precisely, they use the function 

S{z) = 2G{z) + V{z) (4.18) 

rather than G(z) itself. This function S(z) coincides with the 5-function in 
the theory of dispersionless integrable hierarchies 



3 This amounts to the G-function in the large- N limit of random matrices. 
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4.3 Thermodynamic limit for 5D theory 

In the case of the 5D theory, we set the parameters q and Q in an ^-dependent 
form as 

q = e- Rh , Q = (RA ) 2 , (4.19) 

where R and A are positive constants with RA < 1, and rescale s and t^s 

as 

s ->■ h^s, t k ->■ /TV (4.20) 

Assuming the same properties for the rescaled density function pih^u), 
one can derive an explicit form of the thermodynamic energy functional 
Sg^ ^[p^]. Note that the kernel function g 5D in the rescaled coordinate 
behaves as 

g^ih-'u) = tT 2 gfi{u) + 0(0- (4-21) 

gf^{u) is a somewhat complicated function; for our purpose, it is enough to 
know that this function can be characterized by the following conditions: 

Ru/2 _ -Ru/2 

= kg ! — > »S'(0) = 0, »S(0) = 0. (4.22) 

Hence the quadratic part of has the following asymptotic form: 

n poo 

/ / dxdyg 5D (\x - y\)Ap(x)A(y) 

J J —OC 

p poo 

= h 2 dudvgf^u - v\)p^\u)p w '(v) + O^" 1 ). 

J J —oo 

On the other hand, the building blocks of the linear part behave as 

log? 
6 

and 



/OO T~> POO 

dx{x - r 1 s) 3 Ap(a;) = h~ 2 ^ / du(u - s) 3 p^'(u) + C^TT 1 ) 
-oo ^ J —oo 



f 

J — < 



/OO / °° ^ 
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The energy functional thus turns out to be expressed as 

4\h-h[p\ = ^ 2 <? ( °V 0) ] + 0{h-\ (4.23) 

where 



C ( °V 0) ] = / / dudvg^u - v\)pW(u)pW(v) 

J J— oo 

p roo poo I 00 1 — Rku\ 

+ | J du(u - sfp^u) +J dui£ p (0) '(«). (4.24) 



v fc=l 



4.4 Riemann-Hilbert problem for 5D theory 

Assuming the same "one-cut" condition (14.91) as in the case of the 4D theory, 
one can reduce the minimizing problem for S^°\p^} to solving the integral 
equation 

(m <m<M!), (4.25) 

where 

00 

V(u) = Y,t k e- Rku , 

k=l 

under the constraints (14. 3 j) and the support condition 

p(°)'( M ) = for u £ [w , Mi]. (4.26) 

To convert this variational equation to a Riemann-Hilbert problem, we 
introduce the complex analytic function 



W(z) = j dvgS'(z-v)pW(v) 
with the complexified kernel 



R(z-v)/2 _ „-R[z-v)/2 
S 5 D (*"*) = lot, < ^ 



(0)/// < - / e 
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W(z) is a holomorphic function on the z-plane with discontinuity along the 
intervals [uo + 27rin/R, Ui + 27rin/R], n G Z. Moreover, iy(z) + i?z/2 is peri- 
odic with respect to the translation z i— >■ z + 2iri/R. The once-differentiated 
form 

pp £ d^s»(h - ^Dp (0) » + + ^ - o 

(«o<«<«i), (4.27) 
of the variational equation (I4.25P can be converted to the equation 

W(u + iO) + W{u - iO) = -R(u- s) - V'(u) (u < u < m) (4.28) 
for W{u ± iO). This equation is supplemented with the equation 

W (u + iO) - W(u - iO) = i ° ( M>n i)' ( 4 _29) 

I — 2rci [u < uq) 

off the support of p^' . The minimizer p^\u) can be read out from W(z) by 
the same formula as (I4.17p . 

The constraints (14.31) can be translated to boundary conditions of W(z) 
at infinity. Unlike the 4D theory, W(z) exhibits different boundary behaviors 
as Re z — > — oo and Re z — > +oo: 



As Re z — oo and ± Im z > 0, 

W{z) = R ^ Z ~ S > =f m + log(i?A ) + 0{e Rz ). ( 1.30) 
As Re z — >■ +oo, 

W(z) = - " v ~ + log(^Ao) + 0(e- Hz ). (4.31) 



These conditions are derived from (14.31) as follows. As Rez — >■ — oo, the 
kernel function behaves as 

9$?{* ~ «) = ±m- log(i?A ) + 0(e**). 

Consequently, as Re 2; — > —00 under the condition ±Imz > 0, W(z) behaves 
as 

dv + — ± ttz - log(i2Ao)J P (0) » + O(e^) 

=F 7ri + log(.RAo) + 0(e i 



— — — — — — T 77/ ■ 1< IP ( /? V, 1 O'l I '"' ' ' 

2 2 
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In much the same way, as Re z — > +00, 



hence 

W{z) = £ dv(^-^- log(itAo)) p {0) » + 0(e- Rz ) 
= -^ + ^ £ + log( J RAo) + 0(e- fc ). 



Thus W{z) turns out to satisfy (OBI) . (ICTL f TOOj) and flOTjl . We have 
to solve these equations under the condition that W(z) + R(z — s)/2 be a 
periodic function with respect to the translation z H- z + 2ni/R. In other 
words, this is a Riemann-Hilbert problem on the cylinder C/(2iri/R)Z with 
a single cut along [u , 



5 Solution of Riemann-Hilbert problem for 
4D theory 

5.1 Construction of solution when t = 

When t = 0, we can solve the Riemann-Hilbert problem deductively (namely, 
without any heuristic consideration) as follows. This solution is a prototype 
of the solution in the case of t 7^ 0. 

In this case, the equations (14.141) and (I4.15P for W{u ± iO) read 

W{u + iO) + W{u - zO) = for u < u < u x 

and 



W(u + iO) - W(u-zO) 



for it > Mi, 
-27ri for u < Uq. 



One can readily see from these equations that the identity 

e W(u+iO) 1 e -W{u+iO) _ e W(u-iO) _|_ e -W(u-iO) 

holds along the whole real axis. This means that the complex function e w ^ + 
e -W(z) no discontinuity, hence becomes a holomorphic function on the 
whole plane C. Moreover, by (I4.16p . this function has the asymptotic form 

e W{z) +e -W(.z) = Z _^_± + 0{z -l ) 
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as z — > oo. By Liouville's theorem in complex analysis, the 0(z l ) term 
disappears and one obtains the identity 

A 

Setting 

y = e~ w{z \ 

one can rewrite the last identity as 

y + y~ x = ^-- (5.2) 

( 15 .2p may be thought of as the equation of a complex algebraic curve, 
which is exactly the Seiberg-Witten curve of the undeformed 4D U(l) gauge 
theory. One can solve (I5.2p for y = y(z) as 

V(z) = Z ~° P(z) = (z - sf - AAl (5.3) 

Mo and Ui are determined to be the endpoints of the interval of R where 
P(u) < 0, namely, 

u = s- 2A , u x = s + 2A . (5.4) 



P(z) is understood to be the branch on the cut plane C \ [uq, u\) such that 



y/P(z) = z — s + 0{z ) as z — v oo, 

(5.5) 

±Ihia/ P[u ± zO) > for u a < u < U\. 

Thus a (unique) solution of the Riemann-Hilbert problem can be obtained 
explicitly as 

W{z) = -logy(z). (5.6) 

As a cross check, let us examine the behavior of y(z) and W(z) along 
the real axis. y(z) has discontinuity along [mo,mi], and takes real values on 
R \ [uq, Mi]. The boundary values y{u ± zO) for u < u < u\ are unimodular 
(\y(u±iQ)\ = 1), because they are mutually conjugate imaginary solutions of 
( 15.21) . Bearing these facts in mind, one can trace the behavior of y(u ± iO) as 
u decreases along the real axis. When u starts from the right side of ui, y(u) 
decreases towards the value y{u\) = 1. When u passes u\, y{u) bifurcates 
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into the two imaginary numbers y(u ± iO), which move on the unit circle 
in opposite directions towards —1 until u reaches uq. Accordingly, when u 
decreases from u — U\ to u — uq, aigy{u ± iO) varies from to ±7r. For 
u < uq, y(u) is negative and argy(u ± iO) takes the constant values ±7r. 
Translated to the language of W(z), this behavior of a,rg y(u ± iO) implies 
exactly that (14.151) holds. By (14.171) . one can read out an explicit expression 
of p* (u): 

(u > Ui), 
pi 0) (w) = { ±&rgy(u±tO)/ir (u <u< u x ), (5.7) 

1 (u < U ). 

Let us mention that one can further rewrite this result to the celebrated 
"arc-sine law" [SUES]. 

5.2 Construction of solution when t ^ 

When t 7^ 0, we resort to a heuristic method. Namely, we seek a solution of 
the Riemann-Hilbert problem as a deformation of the foregoing solution. A 
naive way will be to add —V'(z)/2 as 

V'(z) 

W(z) = -]ogy(z)--P-, 



which does satisfy the equations (14. 14j) and ( 14. 15j) for W(u±iO), but not the 
boundary condition (14.161) . To fulfill the boundary condition, we modify this 
naive form as 

W(z) = -logy(z) + N(z)y/P(3 - (5.8) 

where 

• y{z) is a solution 



y(z)= Z ^A P{z) = {z-Pf-AA\ (5.9) 



of the equation 



y+y-^tJL (5.10) 



of the deformed Seiberg-Witten curve. 
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• P and A are functions (3(s, t) and A(s, t) of s and t that reduce to the 
previous values at t = 0: 

P(s,0) = s, A(s,0)=A (5.11) 

• Uq and Mi are the zeros of P(z): 

u = /3 -2A, Ul = (3 + 2A. (5.12) 

• a/ P(z) is the branch on the cut plane C \ [uq, Ui] such that 

^pjz) = z-P + 0(z~ asz^oo, 

(5.13) 

±Imy P[u ± iO) > for uq < u < u\. 

• N(z) is a linear combination 

oo 

N(z)=J2 t kN k (z) 
fc=i 

of polynomials Nk(z), k — 1, 2, • • • , in z. 

(I4.14p and (14.15P are already satisfied by the ansatz (15. 8p . To fulfill the 
boundary condition (I4.16p . we choose Nk(z) to satisfy the condition 

k k ' 

N k (z) - -^== = 0(z- 1 ) as^oo. (5.14) 

Nk(z) is uniquely determined by this condition as 

N 1 (z) = 0, N k (z) = ^(z k ~ 2 + Cl z k ~ 3 + ■ ■ ■ + c,„ 2 ) for k>2, (5.15) 
where ci, c 2 , • • • are the coefficients of the expansion 

! z' 1 + Ci^T 2 + c 2 z' 3 H , 



and we set c = 1 for convenience. N k (z) can be also expressed by using a 
contour integral as 

JV*(*) = — </ dX ^ + (5.16) 

kK > 2mJ c x-z2^P[x) l^fpjz) 1 ' 
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where C is a simple closed curve that encircles the interval [uq, Ui] anti- 
clockwise and leaves z outside. Rewriting W(z) as 

00 / U 7 k-1 \ 
W(z) = - logy(z) - J> -j== - N k {z) ^Pjzj, 

one can see that 

W{Z) = - log Z '^ + 2 f^ - £ ^(C^ + C fc Z- + ■ • • ) y/P® 

k=l 

= -log^--^^— + U +0(2; ) 

fc=i V fc=i / 

as 2 — > 00. Matching this expression with (I4.16p . one obtains the equations 

lo A_|,t^ tl = 0| (517) 

/3 - f; Hfc(Cfc ~ = a (5.18) 
?c=i 

for /3 = (3(s, t) and A = A(s, t). 

By the implicit function theorem, these equations do have a solution in 
a neighborhood of t — that satisfies the initial condition (15. lip . Though 
we omit details, ( I5.17P and (I5.18P are equivalent to the equations that are 
derived by Marshakov and Nekrasov by their method based on the S'-function 

(gup Hang. 

Once W(z) is thus determined, one can use (14.171) to obtain an explicit 
expression of p* (u): 

(m > Mi), 
p! 0) (m) = { ±argM(M±zO)/7r + Ar(M)v/|P(M)| (m < u < Mi), (5.19) 

1 (m<m ). 

Note that the structure of (15. 7p is retained except that a new term propor- 
tional to a/|-P(m)| is added. 

5.3 Rewriting solution in terms of Lax function 

By construction, y = y(z) is (a branch of) the inverse function of 

z = z(y) = P + Aiy + y- 1 ). 
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We now interpret z(y) to be a long-wave limit of the well known Lax operator 

£ = a(s)e 9s + b(s) + a(s - l)e~ d % d s = d/ds, (5.20) 

of the Toda lattice [32], where e ±ds are the shift operators, namely, e ±ds f(s) = 
f(s ± 1). Note that one can rewrite the last term of (I5.20p as 

a(s-l)e~ ds =e~ ds -a(s). (5.21) 

In the long-wave (or "dispersionless" ) limit, e 9s turns into a c-number variable 
p (see Section 6 for details). We identify y with this variable p. 

To pursue this analogy further, we introduce the truncation notations 

J2 an y n ) =J2 an y n > [J2 a ^y n ) = J2 an y n > 

\ n / >0 n >0 \ n / <0 n <0 

^ a n y n J = a m (m G Z) 

v. n / m 

for Laurent series of y; the same notations are used for difference operators 
as well. The time evolutions of the dispersionless ID Toda hierarchy are gen- 
erated by the Laurent polynomials (z(y) A: ) > o — (z(y) k ) <0 , k = 1, 2, • • • , which 
are dispersionless limit of the generators (£ fc )>o — (£ k )<o of time evolutions 
of the ID Toda hierarchy |33] . 

We now show that the foregoing polynomials Nk(z) are closely related to 
these Laurent polynomials of y. Note that (z(?/) fc ) >0 — (z(y) k ) < o is a linear 
combination of y l — y~ l , I = 1, 2, • ■ ■ , e.g., 

(z(y)) >0 -(z(y)) <0 = Aiy-y^ 1 ), 
(z(y) 2 ) >0 - (z(y) 2 ) <0 = A% 2 - y- 2 ) + 2A/3(y - y' 1 ), ■ ■ • . 

Since y k — y~ k can be factorized as 

y k - y- k = {y k - 1 + y k ~ 2 + ■■■ + y~ k+2 + y- k+1 )(y - y~ l ) 

and y(z) and its inverse satisfy the relations 



y(z)+y(z)-' = Z - 1 l, y(z)-y(z)- 1 = ^S.. 

i z {y) k )>o — ( z (y) k )<o can be expressed as 

(z( y ) k ) >0 -(z( y ) k ) <0 = Q k (z)^P{i), 
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where Q k (z) is a polynomial in z, e.g., 

(z(y)) >0 -(z(y)) <0 = y/P{z), 
(z(y) 2 ) >0 -(z(y) 2 ) <0 = (z + (3)^P(z), 

The polynomial Q k (z) can be determined as follows. Rewrite the last identity 



as 



z(y) h )>o - (z(y) k )< n ( . 

— Wk{z) 



and note that 

LHS = — ^= + Oiz^ 1 ) as z ->■ oo. 

Truncating the negative powers of z yields the identity 

Q k (z) = z h ~ x + c lZ k - 2 + ■■■ + c fc _i. 
One can thus derive the fundamental formula 



(z(y) k ) >0 - (z(y) k ) <0 = {z k ~ l + Cl z k - 2 + ■■■ + c k . x )y/P{£). (5.22) 
Comparing this formula with (15.1 5p . one can readily see that 

Consequently, the second term on the right hand side of (15.81) turns out to 
be expressed as 

N{z)</P{zj = -f {('to^U " (Avf-'U) • (5-24) 
fe=i 

Moreover, since 

(z(y) k ) >0 ~ (4v) k )<o = z k - {z{y)% - 2(z(y) k ) <0 , 
one can rewrite (I5.22p as 

(z(y) k ) + 2(z(y) k ) <0 = z{yf - (z(y) k ) >0 + (z(y) k ) <0 

= z k - (z^ 1 + Cl z k ~ 2 + ■■■ + C fc _x)/P(i) 

= {c k z~ l + c k+1 z- 2 + ■ • -)y/P(zj 
= c k + (c fc+ i - (3c k )z~ l + 0(z~ 2 ). 
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Since z{y)~ l = (Ay) -1 + 0(y~ 2 ) as y — >■ oo, one can re-expanded the right 
hand side in powers of y~ l and pick out the coefficients of y° and y~ x . This 
leads to the identities 

(z(y) k ) = c k , 2(z(y) k U = Ck+1 -^ Ck . (5.25) 

One can thereby rewrite (15.171) and (15.181) as 

logf -£^(^) fe - 1 )o = 0, (5.26) 
A ° k=i 2 

oo 

/^-A^H^^- 1 )^*- (5.27) 
fc=i 

In the next section, we derive the same equations from generalized string 
equations. 



6 Generalized string equations for 4D theory 
6.1 2D Toda hierarchy 

Let us briefly recall the construction of the Toda hierarchy |1 7] @ (referred to 
as the "2D" Toda hierarchy to distinguish it from the more classical "ID" 
hierarchy[33j). 

This integrable hierarchy is formulated in terms of two Lax operators L 
and L, which are difference (or, so to speak, "pseudo-difference") operators 
with respect to the lattice coordinate s. We now choose the "symmetric 
gauge" : 

oo oo 

L = ae 9 * + u n e {1 - n)d ° , L = a~ l e 9s + ^ uj l+n ^ , 

n=l n=l 

where a, w„'s and u„'s H are dynamical variables that depend on s and the 
two sets of time variables T = (Ti, T2, ■ • • ) and T = (7\, T2, • • ■ ). Let us 

4 Since the first paper [T7] was published, the notations for the Toda hierarchy have 
been considerably changed. Our notations are mostly based on the review [18J . 
5 The bar does not mean complex conjugate. 
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introduce the truncation notations 



\ n / >0 n >0 \ n / <0 "<0 

J2 a ne n9s ) =a m (meZ). 



7)1 



for difference operators. 

The time evolutions of the Lax operators are generated by Lax equations 
of the form 



dL 

m 


= [Bk- 


L], 


dL 

df k 


= [Bk- 


L], 


dL 

m 


= [Bk- 


L], 


dL 
df k 


= [Bk, 


L], 



(6.1) 



where 

B k = (L k ) >0 + i(L fc ) , B k = (L~ k ) <0 + \{L-%- 
These Lax equations are supplemented by another set of Lax equations 

^ = [S*,M], 
Oik oik /„ s 

dM 8M - - 1 ' 

jgr = [*,*], Wt =M 

and the twisted canonical commutation relations 

[L, M] = L, [L, M\ = L (6.3) 
for the Orlov-Schulman operators 

oo oo oo oo 

M = J2 kT kL k + S + v nL~ n , M = -J2 k TkL~ k + S + v n L n . 

k=l n=l fc=l k=l 

6.2 Dispersionless 2D Toda hierarchy 

The dispersionless Toda hierarchy[15j is a long-wave limit of the Toda hier- 
archy. In this limit, as briefly mentioned in Section 5, the shift operator e " 
is replaced by a momentum-like variable p. p and s then become a (twisted) 
canonical coordinate system of a 2D phase space of classical mechanics. This 
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procedure is justified by applying the idea of semi-classical (or quasi-classical) 
approximation in quantum mechanics to the Toda hierarchy [T5] . 

To derive the dispersionless Toda hierarchy, one starts from an ^.-dependent 
formulation of the Toda hierarchy, in which e ds , d/dTk and d/dTk are re- 
placed by e hds , hd/dTk and hd/dTk- The Lax and Orlov-Schulman operators 
are assumed to take an ^-dependent form as 



ae 



oo oo 

hds 



+ ]T uj 1 -^ , L = a-V^ + «ne (1+n)Ms , 

n=l n=l 

oo oo oo oo 

m = j2 kT kL k + v ^ L ~ n ' M = -J2 kfkL ~ k + s + Yl 

k=l n=l k=l k=l 

where a, u n ,u n and v n ,v n depend on h as well as (s,T,T) and have semi- 
classical expansions 

a = a (0) + ha {1) + ■■■ , 
u n = + hu^ + ■■■ , u n = u < £ ) + fifiW + • • • , 

Wn = V® + + • • • > = ^ + ^i X) + • • • 

as ft — ?• 0. They obey the Lax equations 

H— = [B k ,L\, ft_ = [B t ,L], 

(equations of the same form for M) • • • , 
and the twisted canonical commutation relations 

[L, M] = HL, [L, M] = hi. (6.5) 
One can derive, from these equations, the Lax equations 

(equations of the same form for M.) • • • , 
and the twisted canonical relations 

{£,M} = £, {£,M} = £ (6.7) 
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for the Laurent series 



£ = a^p + Y C = a^-V + £ <V + ", 

n=l n=l 

oo oo oo oo 

M = J2 kT k£ k + S + V n )jC ~ n > M = ~Y kf kC~ k + S + J2 VTC 1 

k=l n=l k=l k=\ 

of the new variable p. B k and B k are determined by £ and £ as 

B k = (£ k ) >0 + \{£ k \, B k = (£- h ) <0 + l(^ k )o, 

where ( ) >0 , ( )<o and ( ) are now understood to be truncation opera- 
tions for Laurent series of p: 

Y a ^p n ) = Y anpn > ( Y anpH ) = Y anpU > 

\ n / >0 ™>0 \ n / <0 ra <0 

Y a nP n I = a m (me Z). 
\ n / m 

{ , } denotes a Poisson bracket specified below. These equations are fun- 
damental building blocks of the dispersionless Toda hierarchy. 

This procedure may be thought of as a kind of "classical limit" from 
quantum mechanics to classical mechanics. The Laurent series £,M,£, M, 
are classical counterparts (referred to as the Lax and Orlov-Schulman func- 
tions) of L, M, L, M, in which e 8 turns into the c-number variable p, and 
the coefficients are replaced by the leading terms of the ^-expansion. The 
Poisson bracket { , } is defined as 

{F G \ = P (<!L!K-<!L<*±\ 

\dp ds ds dp ) 

for functions on the 2D phase space (p,s). In particular, p and s obey the 
twisted canonical relation 

{p,s}=p, (6.8) 

which is exactly the classical limit of the quantum mechanical commutation 
relation 

[e hd °,s] = he hds . (6.9) 
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As a technical remark, let us note that one can derive the ^-dependent 
Toda hierarchy from the ^-independent formulation by rescaling the variables 
s,T,t as 

s-tfrh, T k ^H~ l T k , f k ^fr x f k (6.10) 
and simultaneously rescaling M and M as[§ 

M -> fT x M, M -> n~ x M. (6.11) 



This remark is very important for our purpose. Note that (I6.10P amounts 
to (14. ip and (14.201) in the prescription of thermodynamic limit, [j By the way, 
the energy functional (13.181) and (13.191) depend on h (through the relation 
q = e~ Rh in the 5D case) before the rescaling. This is not a contradiction, but 
necessary to obtain a non-trivial thermodynamics limit. In a similar sense, 
to obtain a non-trivial solution of the dispersionless Toda hierarchy, one has 
to start from a h- dependent solution of the ^-independent Toda hierarchy, 
and let H — >■ after rescaling the variables as (16.101) and (16.111) . Moreover, 
for the rescaled solution to have a limit as H — > 0, one has to choose the 
^-dependent solution carefully . 

6.3 Generalized string equations 

Since we shall not return to the Toda hierarchy in the rest of this paper, let 
us use the same notations L, M, L, M, rather than C, J\4, C, Ai, for the Lax 
and Orlov-Schulman functions of the dispersionless Toda hierarchy. They 
are assumed to be Laurent series of the form 

oo oo 

L = ap + 2^ u nV X ~ n i L = or x p + 2J u„p 1+n , 

n=l n=l 

00 00 00 00 

M = kT kL k + s + v nL~ n , M = -J2 kf kl- k + s + Y, v n l n . 

k=l n=l k=l k=l 

More precisely, L and M are understood to be Laurent series in a neighbor- 
hood of p = 00, and the other two in a neighborhood of p = 0. To consider 



6 Only after this rescaling, the twisted canonical commutation relations of the Lax and 
Orlov-Schulman operators take the correct form (|6.5p . Actually, this ad hoc reasoning can 
be justified more rigorously in the language of an auxiliary linear problem [18] . 

7 (|4.ip looks apparently different from (|6.11l) . but actually, they are substantially the 
same. A careful inspection shows that and $fc(/i, s) are defined to have an excess factor 
of h. Because of this factor, tk's have to be rescaled by h~ 2 rather than H . 
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generalized string equations, these four Laurent series are further assumed 
to have a common domain of convergence, say, r < \p\ < r^. @ 

Solutions of the dispersionless Toda hierarchy in a general position [18] 
are characterized by functional equations of the form 

L = f{L, M), M = g(L, M), (6.12) 

where / = f(z,w) and g = g(z,w) are arbitrary functions of two variables 
that satisfy the symplectic condition 

df Adq dz A dw 

J y = • (6.13) 

/ z 

If the functional equations ( I6.12p have a solution, namely a quartet (L, M, L, M) 
of Laurent series of the form shown above, they automatically give a solu- 
tion of the dispersionless 2D Toda hierarchy [18]. Because of their origin in 
string theories [191 EDJ [211 1221 [23] , (16. 1 2[) are referred to as "generalized string 
equations" . 

The goal of our consideration to the end of this section is to identify 
the solution ( 15. 8p of the Riemann-Hilbert problem with a solution of the 
particular generalized string equations 

L - l' 1 

L = L-\ L~ l M — log — = —LM + log , (6.14) 

ao a 

where ao is a positive constant that will be eventually identified with A . Let 
us mention that we have been unable to derive (16.141) from the first principle, 
namely, from properties of the partition function. We reached ( I6.14p by guess 
work seeking analogy with the 5D counterparts. In the case of the 5D theory, 
we know a derivation from the first principle (see Section 7). 

It might be, however, possible to derive (I6.14p from a suitable modifi- 
cation of the Eguchi-Yang model [3"6] . This random matrix model has a 
logarithmic term of the form Tr(XlogX — X) in the potential, which will 
lead to the logarithmic terms of ( I6.14p . See also the work of Nekrasov and 
Marshakov [12, 13J on the origin of logarithmic terms in their approach to 
the thermodynamic limit. 

A few technical remarks on (I6.14p will be in order. 



8 One can also consider a genuinely algebraic framework in which the coefficients of the 
Lax and Orlov-Schulman functions are meaningful as formal Laurent series of p. In such 
a formulation, the coefficients u n ,v n ,u n , v n are treated as formal power series of T and 
T, whose coefficients are functions of s in a suitable class. 
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• (I6.14p do not take the form shown in (16.121) . This is not a serious 
problem. What is crucial is not this standard form itself but the sym- 
plectic condition (I6.13p . This condition ensures that the map (z, w) h-> 
(f(z,w), g(z,w)) connecting (L,M) and (L,M) is symplectic with re- 
spect to the symplectic form dz/z A dw. In the case of (16.141) . the 
relevant map (z, w) h-> (z, w) is defined implicitly by the equations 

Z Z 

z = z~ x , z~ x w — log — = — zw + log . (6.15) 

ao ao 

It is easy to show that this map is indeed symplectic: 

dz A dw dzAdw 



• The logarithmic terms in (16. 14j) have multi-valuedness in the vicinity of 
p = and p — oo in which the Lax and Orlov-Schulman functions are 
assumed to have Laurent expansion with respect to p. This problem 
can be settled by a simple trick that we shall present in the course of 
solving (I6.14p . 



6.4 Solution of generalized string equations 

Let us construct the solution of (16.141) that is to be identified with the solution 
of the Riemann-Hilbert problem. Apart from the trick for the logarithmic 
terms (which resembles a trick in the Eguchi-Yang model [22]), one can P ro " 
ceed in much the same way as the case of generalized string equations in the 
large- N Hermitian random matrix model [37J. 
The first equation of (16.141) implies that 

£ : = L = L- 1 



is a Laurent polynomial of the form 

£ = ap + b + ap' 1 , (6.17) 

where a and b are functions of s, T and T. This is nothing but the Lax 
function of the dispersionless ID Toda hierarchy, namely, the long-wave limit 

of (Km . 

As regards the second equation of (I6.14p . we add logp to both hand sides 
and reorganize this equation as 

ajt:=L -1 M-log^— = -ZM + log— ^. (6.18) 
a a 
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The modified logarithmic terms, unlike log(L/a ) and log(L 1 /a ), can be 
expanded to Laurent series of p as 

Lp~ x a ( b _i _ 2 

log = log h log 1 + ~p + p 

ao Oq V a 



a b , / 
log — + -p" 1 + 1 
a a \ 



2^ J 



and 



, L p . a , { , b 
log = log h log 1 + -p + p x 

ao ao \ a 

, a 6 A & M 2 

= log h ~p + 1 - — p + 

a a \ ia^y 

In particular, 59T itself is a Laurent series of p. 
(I6.18P splits into following two equations: 



Tt = IT X M — log — — , (6.19) 
a 

SOI = -LM + log^-^. (6.20) 
a 



Substituting 



fc=l n=l 

one can rewrite (I6.19P as 

oo oo , 

Tl = V fcTfcL^" 1 + sL- 1 + V v n L- n - 1 - log — - -p" 1 + • • • . 
^— ' ^— ' a a 

fc=l n=l u 

The ( ) >Q part and the p° part give the equations 

oo oo 

(Tl) >0 = kT k (L k ~ 1 ) >0 = miZ^U, (6.21) 

k=2 k=2 

oo oo 

(Wl) = V kT k {L k ~\ - log - = V kT k {Z k ~\ - log -. (6.22) 

fc=l 

In the same way, f l6.20p reads 



oo oo , 

2JT = V A;f fc Z fe+1 - sL - y v n L n+1 + log - + -p + 



35 



and the ( ) <0 part and the p° part give the equations 

oo oo 

(9Jt) <0 = kf k {L l - k ) <0 = J2 kf^- 1 )^, (6.23) 



k=2 k=2 

oo 



{Wt) = J2 kUL^o + log - = J2 kfki^h + log -. (6.24) 

Equating the two expressions (16.221) and (I6.24p of (2Jl)o, one obtains the 
explicit expression 

Wo = E /C(Tfc 2 +Tfc) (£ fc - 1 )o (6.25) 
fe=i 

of (971) o along with the equation 



oo 



log- = E^P Z (^ 1 V (6-26) 
From flOTj) . (E23D and ( 1Q5j) . one finds that 



OO OO OO 7 /rrn , rxn \ 

m = £ kT k (2 k -% + £ feT*^" 1 )^ + E c , V 'V ( 6 - 27 ) 



Lastly, the ( ) <0 part of (I6.19p and the ( ) >0 part of (I6.20p yield some 
more equations. The equations from the p~ x terms of (I6.19P and those from 
the p terms of (I6.20p give the same equations 



, oo 
S — 



+ £ k(T k - T,.)^- 1 )-! = 0. (6.28) 



k=2 



This equation supplements (16.261) to determine a and b as functions of (s, T, T). 
The equations from higher orders of p ±Y read 

^— + (linear combination of V\, ■ ■ ■ , t> n _i) 



a n+l 



= -Y J KT k -T k ){Z k - l )^ l +{\og(Lp- 1 /a,))_ n _ l , (6.29) 

k=l 

and 



a n+l 



+ (linear combination of V\, • ■ ■ , v i) 

oo 

= -Y^k(T k -T k )(2 k - 1 ) n+1 + (log(Z-V«o)) n+1 . (6.30) 



fc=i 
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These equations determine v n and v n recursively once a and b are obtained. 
One can thus construct a (unique) solution of the generalized string equations 

dHUD. 

6.5 Identification 

As remarked above, the first equation of (16.1 4p is a reduction condition to 
the dispersionless ID Toda hierarchy. If this equation holds, then 

B k + B k = £ fc , (6.31) 

hence 

H + §- = = 0. (6.32) 

This implies that a and b (and, actually, all other dynamical variables) de- 
pend on T and T through the difference T — T only: 

a = a(s,T -T), b = b(s, T — T). (6.33) 

This explains why T k and show up in (16.261) and (16.281) in the form of the 
difference T k —T k . 

Let us now restrict the time evolutions to the anti-diagonal subspace 

T fc = -f fc = |, fc = l,2,-... (6.34) 

The reduced Lax function £ therein satisfies the standard Lax equations 
r}£ 1 1 

— = {A k} £}, A fc = -(£*) >0 - -(£ fc )<o (6.35) 

of the dispersionless ID Toda hierarchy. Moreover, (I6.26p . (I6.28P and (I6.27P 
take the reduced form 

log ^ = ^^ (£fe -i )o) (06) 
a ° k=i 2 

, oo 

^ + ^H fc (£ fc - 1 )_ 1 = 0, (6.37) 

k=2 

^=E^((^ 1 )>0"(^ 1 )<o)- (6-38) 

k=2 
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For (I6.36P and (I6.37P to coincide with (I5.17P and (I5.18p . it is sufficient that the 
building blocks of the Riemann-Hilbert problem and the generalized string 
equations are identified as 

a = A , a = A, b = (3, y = p, z(y) = £. (6.39) 

Comparing (16.381) with (I5.24p . one finds that that W(z) and 0JT are related 

as 

V'(z) 

W(z) + log y(z) + = N{z) y/P(z) = Wl. (6.40) 

Thus (15. 8p can be identified with a solution of the dispersionless Toda 
hierarchy that satisfies the generalized string equations (16.141) . Let us stress 
that the generalized string equations are a nonlinear analogue of (so to speak, 
modern) matrix Riemann-Hilbert problems that are widely used in the theory 
of integrable systems [18]. We thus have a nonlinear (ultra-modern) Riemann- 
Hilbert problem for L, M, L, M on the p-plane. On the other hand, we have 
another (classic) Riemann-Hilbert problem for W(z) on the z-plane. This is a 
very interesting interplay of two Riemann-Hilbert problems of quite different 
natures. 



7 Solution of Riemann-Hilbert problem for 
5D theory 

7.1 Construction of solution when t = 

As in the case of the 4D theory, we can deductively solve the Riemann-Hilbert 
problem in the case where t = 0. 

When t = 0, one can rewrite 04.281) and 04.291) as 

W(u + i0) + R ( U ~ s) ^ + (w(u-iQ) + R ( U ~ s) ^j =0 ioiu <u< Ul 

and 

R(u-s)\ /„. , R(u-s)\ |0 (u>ui), 
2 J \ 2 J I — 2m (u<u ). 

One can see from these equations that the complex function e w ( z )+ R ( z ~ s )/ 2 _|_ 
e -w(z)-R(z-s)/2 ^g^g no discontinuity along the real axis, hence becomes a 
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holomorphic function on the whole z-plane C. Moreover, by the other condi- 
tions in the Riemann-Hilbert problem, this function is periodic with respect 
to the translation + 2iri/R and behaves as 

e W(z)+R(z-s)/2 + e -W{z)-R{z-s)/2 _ e ^ + q/jx 

as Re z — > — oo and 

e W(z) + R(z-s)/2 + e -W{z)-R{z-s)/2 = J_ + RAq + (e~ Rz ) 

RA 

as Re z —7- +00. Therefore, by Liouville's theorem, one can deduce that the 
identity 

e W(z)+R(z-s)/2 + e -W(z)-R(z-s)/2 _ + R ^ _ e ^ &) /y-g 

holds. 

Let us introduce 

y = e -w(z)-R(zs)/2^ f 3 = {l + (RA f) e -^ : A = A oe -^ 
to rewrite the last identity as 

_ -Rz 

Viewed as an equation of a complex analytic curve, f)7.2p defines substantially 
the same Seiberg-Witten curve as previously considered for the undeformed 
5D U(l) gauge theory [14|. 

One can solve this equation for y as 

y = y^) = P ~ e ' R ' 2 R/ m ^ P(z) = (e- Rz -P) 2 -1(RA) 2 . (7.3) 

Mo and u\ are determined to be the endpoint of the interval of R where 
P{u) < 0, namely, 

log(/3 + 2 J RA) 21og(l + i2A ) 



K K , 7 4 n 

log(/3-2 J RA 21og(l-i?A ) 1 ; 

Ul = R = S R ' 
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y/P(z) is the branch on the z-plane cut along the intervals [u + 2irin/ R, ui + 
27rin/R, n e Z, such that 



'P(z) = e~ Rz -0 + 0(e Rz ) as Re z -»■ -oo, 
v/p(i) = -v//3 2 - 4(i?A) 2 + 0{e- Rz ) asRez^+oo, (7.5) 



^plmy P(u ± iO) > for uo < u < Wi- 

Solving 

e -W( Z )-R( Z -s)/2 = y{z) 

for W(z) yields a (unique) solution of the Riemann-Hilbert problem explicitly 

as 

W(Z) = - R{Z ~ S) -logy(s). (7.6) 

As explained in the case of the 4D theory, it is easy to reconfirm directly 
that this function fulfills all requirements of the Riemann-Hilbert problems. 
The minimizer p^\u) has the same expression as the 4D case (15. 7p . 

7.2 Construction of solution when t ^ 

Let us now consider the case where t ^ 0. We seek a solution of the Riemann- 
Hilbert problem by modifying the solution for t ^ as 

W(z) = ~{z -s)- logy(z) + N(z)^Pjz)- (7.7) 

where 

• y(z) is a solution 



y(*) = - e 't„ A V ^ , P{z) = (e" & - (3) 2 - 4(i?A) 2 , (7.8) 



2RA 

of the equation 



a _ -Rz 



of the deformed Seiberg-Witten curve. 

(3 and A are functions (3(s, t) and A(s, t) of s and t that reduces to the 
previous values at t = 0: 

/3(s, 0) = (1 + (it:Ao) 2 )e-^, A(s, 0) = A e- fo . (7.10) 
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u and Mi are the real zeros of P(z): 

log((3 + 2RA) \og((3-2RA) 
u = ^ , ui = ^ • (7.11) 

^JP(z) is the branch on the z-plane cut along the intervals [u + 
2Tcin/R,Ui + 2mn/R], n e Z, such that 



v/p(i) = e~ Rz -(3 + 0{e Rz ) as Re z -> -oo, 

= - - 4(#A) 2 + 0(e~ & ) asRez^+oo, (7.12) 



=F Im y P(« ± zO) > for uq < u < u\. 
• N(z) is a linear combination 

oo 

N(z) =Y^hN k (z) 

k=l 

of polynomials N k (z) in e~ Rz . 

f)4.28p and (14.291) are already satisfied by the ansatz (17. 7J) . Moreover, 
since -P(z), A^(z) and V(z) are periodic with respect to the translation z i-» 
2 + 2m/ R, the periodicity condition of + R(z — s)/2 is also satisfied. 

Therefore it is sufficient to fulfill the boundary conditions (I4.30P and (14.311) . 

To this end, we choose N k (z) to satisfy the condition 

W p -Rkz 

N k {z) + ^-^= = 0{e Rz ) asRe^-oo. (7.13) 



N k (z) is uniquely determined by this condition as 

N k (z) = -?±(e~ R ^ z + Cl e~ R ^ z + ■■■ + cO, (7.14) 
where c±, C2, ■ ■ ■ are the coefficients in the expansion 

' e Rz + Cl e 2Rz + c 2 e 3Rz + ■■■ , 



and we set c = 1 for convenience. Consequently, as Rez — > — oo in the 
domain ±lmz > 0, W(z) behaves as 



R(z - s) 13 - e~ Rz - JP{z) 

W[z) = — — - - log — 



2RA 

k ( Rz i 

^ =- ( (■>.< 4- r /. 



+ E ^ i c ^ RZ + ^ + ie R(fc+1)z + • • • ) VW) 

k=l 



fc=l 
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One thus finds the equation 

^ , A RktkCk „ 

ifc + log — + — f 1 = ( 7 - 15 ) 
A ° k =i 2 

for the first boundary condition (14.301) to be satisfied. 

On the other hand, in the opposite end (Re z — » +oo) of the z-plane, 
N(z), P(z) and y(z) converge to finite values as 

lim N(z) = - y RktkCk ~ 1 



k=l 



lim ^/P(z) = -A(RA) 2 , 

Re z— >+oo 



y (3+ v//3 2 -4(i?A) 2 

yoo := hm y(z) - 



Rez-^+oo 2RA 

and V(z) disappears. Therefore W(z) behaves as 

R( z ~ s ) i , Rkt^k 

k=i 

One thus obtains the equation 



- logyoo + E^T^V^-^A) 2 = log(i?A ) (7.16) 
fc=i 

for the second boundary condition (14.311) to be satisfied. 

Thus the problem reduces to solving (I7.15P and (I7.16P for f3 = (3(s,t) 
and A = A(s, t). The implicit function theorem ensures the existence of a 
solution in a neighborhood of t = that satisfy the initial conditions (I7.10p . 
Note that f T7TT5j) and ( 17715]) reduce to 



p , , A P+y/P- 4(i?A) 2 1 

its + log — = 0, 



A 2i?A RA 

when £ = 0. It is easy to see that (I7.10P give a solution of these equations. 
Lastly, the the minimizer p^\u) turns out to have the same expression as 
the 4D case <^Wj. 
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7.3 Rewriting solution in terms of Lax function 

Let us rewrite the foregoing solution of the Riemann-Hilbert problem in terms 
of the new variables 



Z = e 



-Rz 



Note that the cylinder C/(2iii/R)Z is thereby mapped to the punctured 
plane C* = C \ {0}. P(z) and N k (z) thereby turn into polynomials in Z as 

P(z) = {Z- /3) 2 - 4(i?A) 2 , 
N k (z) = ~^(Z k ^ + Cl Z k ~* + ■■■ + c fc _0. (7 ' 17) 
c„'s may be thought of as the coefficients of the Laurent expansion 
' Z- 1 + c 2 Z- 2 + c 3 Z- 3 + ■■■ 



V / (^-/?) 2 -4(i?A) 2 
at Z = oo, and A^(z)'s are uniquely determined by the condition 

F?k7 k 

N k (z) + — _ = OiZ- 1 ) asZ^oo. 

2 v /(Z-^) 2 -4(i2A) 2 

Another expression of N k (z)'s is the contour integral formula 

. . 1 f dX RkX k RkZ k 

N k {z) 



2ni JcX-Z 2^{X - f3) 2 - 4{RA) 2 2 y/(Z - f3) 2 - 4{RA) 2 

(7.18) 

where C is a simple closed curve that encircles the interval [e~ Sui , e~ Ruo ] 
anti-clockwise and leaves e~ Rz outside. 

The function y(z) can be redefined as a function of Z as 

which is an inverse function of 

Z(y) = f3-RA(y + y~ 1 ). 

This function Z(y) plays the role of a Lax function. |^| As in the case of the 
4D theory, one can derive the following identities: 



(Z(y) k ) >0 - (Z(y) k ) <0 = {Z k - 1 + Cl Z k ~ 2 + ■■■ + c k ^)^P{z), (7.20) 
(Z(y)% = c k , (Z(y) k U = - Ck+1 2 -P Ck . (7.21) 



9 Strictly speaking, it is —Z(y) rather than Z(y) that correspond to the Lax function 
£ (see section 8). 
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One can thereby rewrite the third term of (17.71) as 



N{z)y/p{£) = - £ i^y^ ~ ( z (y) k )<o) (7-22) 



fc=i 

and (17.151) as 



R* + log± + £^(^) fc )o = 0. (7.23) 

fe=i 

Unlike (jOSjl . ( jOSD and fl4T30|) . it seems difficult to interpret fl7~T6|) cor- 
rectly within this framework. (17. 16p is derived to fulfill the boundary con- 
dition (14.31 j) as Rez — > +oo, in other words, Z — > (rather than Z — > oo). 
The situation of the 5D theory thus turns out to be considerably different 
from the 4D theory. 



8 Solution of dispersionless Toda hierarchy 
for 5D theory 

8.1 Generalized string equations 

The goal of this section is two-fold. Firstly, following the method of Section 6, 
we show that the solution f lT.Tf) of the Riemann-Hilbert problem corresponds 
to a solution of the particular generalized string equations 

L = L~\ e~ RM = Q- x L- 2 e m . (8.1) 

Secondly, we derive these equations from hidden symmetries of a tau function 
r s (T,T) of the Toda hierarchy [SJ As we briefly review therein, this tau 
function reduces to the partition function Z^ D (t) when the time variables 
(T, T) are restricted to the anti-diagonal subspace. 

Let us remark that (18.11) can be converted to a form that resembles 
the generalized string equations (I6.14p of the 4D theory. Rewrite the term 
Q~ l L~ 2 on the right hand side of the second equation as 

-17-2 L L~~ 



Q~ L L 



Q 1 / 2 Q 1 ' 2 ' 



move the factor L/Q 1 ^ 2 to the left hand side and take the logarithm of both 
hand sides. One thus obtains the equations 

L = L~\ -R M -\og-^ = RM + \og^ (8.2) 
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with remarkable similarity with (I6.14p . 

Unfortunately, (I8.2p turn out to be not strong enough to characterize the 
solution of the Riemann-Hilbert problem uniquely. As one finds in the course 
of solving these string equations, (I4.28p . (I4.29P and (14.301) are correctly en- 
coded in the generalized string equations, but the second boundary condition 
f!4.3ip is missing therein. Because of this, our solution of the generalized sting 
equations contains arbitrariness. This arbitrariness can be fixed by imposing 
f)4.3ip by hand, but its status in the framework of the dispersionless Toda 
hierarchy is obscure. 

Thus, although looking very similar, the 5D generalized string equations 
( 18.21) are drastically different from the 4D version (I6.14p . We can point out 
technical reasons leading to this difference, but a true reason is beyond our 
scope. 



8.2 Solving generalized string equations 

The method for solving (18. 2p is parallel to the case of the 4D theory. 
The first equation of (18. 2p again leads to the Lax function 

£ := l = L- 1 = ap + b + ap- 1 (8.3) 

of the dispersionless ID Toda hierarchy, a = a(s,T,T) and b = b(s,T,T) 
are arbitrary at this stage. We now derive equations for these functions from 
the the second equation of (18.21) . 

Adding logp to both hand sides, one can rewrite the second equation of 
(IQ1 as 

m .= _ RM - log ^ = RM + log , (8.4) 
and separate this this equation to two equations: 

9JI = —RM — log -^p/jf) (8.5) 
m = RM + log (8.6) 
The ( ) >0 parts of (I8.5P and the ( ) <0 part of (I8.6P give 

oo oo 

(m) >0 = -J2 RkT k (2 k ) >0 , (M) <0 = -J2 RkT k (£ k ) <0 . (8.7) 

k=l k=l 
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The ( ) parts of (I8.5P and (I8.6P read 

CO 

(9rt) = -Rs - log-^ -^^(£*) , (8.8) 
^ k=i 

oo 

(OJl)o = Rs + log - J] Rkf k (£ k ) . (8.9) 
^ fc=i 

Equating these two expressions yields the explicit expression 

Wo = _^ ^(n + T fc ) (£fc)o (gio) 
fc=i 

of (9Jt)o and the equation 

Rs + lo g ^ + ± Rkin - Tt) (£*) = . (8.11) 



k=l 



Note that (18. lip amounts to (I6.26P in the case of 4D theory. Thus one finds 
that 

Tl = - £ RkT k (2 k ) >0 - jr RkT k (Z k ) <0 - £ Rk{Tk 2 +Tk) (Z k ) . (8.12) 

fc=l k=l k=l 

On the other hand, from the ( ) <0 part of f)8.5p and the ( ) >0 part of 
(18. 6p . one obtains the recursive formulae 

— - + (linear combination of v%, • ■ ■ , t> n _i) 

oo 

= -Y,Rk(T k -T k )(Z k )_ n - (logOLp-VQ 1 / 2 ))^ (8.13) 

k=l 

and 

+ (linear combination of Vi, • ■ • ,v±) 



- RKTk ~ f k ) (£ fc ) n - (log(Z-W /2 )) n (8.14) 



k=l 



of v n and v n for n = 1, 2, • • • . Note that these equations are all that one can 
derive from the ( )<o part of (I8.5P and the ( )>o part of (I8.6p . There is no 
equation that amounts to (I6.28P in the case of the 4D theory. 



46 



Thus, unlike the case of the 4D theory, we have only one equation ( 18. lip 
for the two unknown functions a and b. Therefore these functions remain 
to be fully determined. This is the aforementioned "arbitrariness". In this 
sense, these generalized string equations ( 18. 2 p are incomplete. 

From a genuinely technical point of view, this phenomenon stems from 
a delicate difference of the structures of (I6.14p and ( 18. 2p . In the second 
equation of ( I6.14p . M and M are multiplied by L^ 1 and L. Because of 
the presence of these multipliers, when this equation is split into (I6.19P and 
(I6.20p and expanded in powers of p, the p~ x terms in (I6.19P and the p terms 
in ( I6.20p yield exceptional equations. These equations, which are mutually 
equivalent, are nothing but the second equation (16. 28j) for a and b. In the 
case of the second equation of (I8.2p . there is no such exceptional terms, hence 
no counterpart of ( I6.28p . 

Let us now restrict the time evolutions to the anti-diagonal subspace. As 
the main result of our previous work [8] predicts, the reduced time variables 
tk should be twisted by the signature factor (— l) k as 

T k = -f k = (-l) fc |, fc = l,2,.... (8.15) 
(18.1 ip and (18.121) then become the following equations: 

Rs + log — - 2 + -(^)o = 0, (8.16) 

^ fe=i 

m = _J- tD^Mh ( (£ * )>0 _ ( ^ )<0 ) , . (s.iT) 
fc=l 

(I8.16P coincides with ( I7.15p . if the building blocks of the Riemann-Hilbert 
problem and the generalized string equations are identified as 

Q = (RA ) 2 , a = RA, b = -/3, y = p, Z = e~ Rz = -£. (8.18) 

Note that the sign factor (— l) k in (I8.15P is correlated with the negative sign 
in the last relation of (I8.18P between Z and £. Comparing (I8.17P with (I7.22p , 
one finds that W(z) and DJl are related as 

W{z) + j(z -s) + logy(z) + = N(z)y/P{z) = -M. (8.19) 

These are all that one can derive from the generalized string equations; (I7.16P 
is missing here, and has to be added by hand for a and b to be determined. 
Thus the aforementioned "arbitrariness" is not resolved even in the diagonal 
subspace (18.151) . 
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8.3 Deriving generalized string equations 

We now turn to the second subject of this section, namely, the origin of the 
generalized string equations in hidden symmetries of a tau function of the 
Toda hierarchy. 

As observed in our previous work [8], the partition function of the 5D 
theory is connected with a tau function t s (T,T) of the Toda hierarchy as 

^ 5D (*)=exp (ErrV) 1- s{s+ms+m 

\k=l q J 

XT s (T,T)\ Tk= _ fk={ _ 1)Hk/2(k > iy (8.20) 
This tau function t s (T,T) has the fermionic expression 

t s (T,T) = (s\exp (^T k J k J ^ ex P f-£^ J -fcj I s )' ( 8 ' 21 ) 

where g is the GL(oo) element 

g = q w »/ 2 G_G + Q Lo G_G + q Wo/2 . (8.22) 

It is also pointed out in the same work (loc. cit.) that g satisfies the 
algebraic relations 

J k g = gJ. k , fc = l,2,-... (8.23) 
These relations can be converted to the equations 
d d \ 

+ -=r ) r s (T,T) = 0, fc = l,2,-.., (8.24) 



,dT k dT k/ 

for the tau function, which thereby reduces to a function of T — T: 

t s (T,T)=tI d (T-T). (8.25) 

The function r s 1D (T) is a tau function of the ID Toda hierarchy. 

In the language of the Lax operators [2H [23], (I8.23j) amount to the rela- 
tions 

L k = L~\ A; = 1, 2, • • • , (8.26) 
which reduce to the single equation 

£ := L = L-\ (8.27) 
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This implies that £ becomes a difference operator of the form (I5.20p , namely 
the Lax operator of the ID Toda lattice. One can thus reconfirm that the 
foregoing tau function is a special solution of the ID Toda hierarchy. 

Actually, (I8.23p are merely a small part of a large set of relations satisfied 
by g. Those relations are derived from the following two types of "shift 
symmetries" [S] among the elements 

nGZ 

1 — q K 

of the quantum torus algebra: 

• First shift symmetry 

G_G + V r ^ { G_G + )- 1 = (-l) k vi%. (8.28) 

• Second shift symmetry 

q w °l 2 V^q- w °' 2 = Vi k ~^ (8.29) 

equivalently, 

q Wo/2y(k) q -Wo/2 = y(k-m)_ (g>30) 

(k) 

Using these symmetries repeatedly, one can rewrite Vm g as 

V^g = Vi k k Wo/2 G-G + Q L »G-G + q w °/ 2 

= q Wo l 2 v£ + ^G_G + Q L »G_G + q w °l 2 

= (-l) k +™qWo/2 G _ G+ y(k+™)QL» G _ G+e W /2 

= (-i) k +™Q- k q w ^G_G + Q Lo V^ +m) G_G + q w °/ 2 

= (-l) fe + m (-l) fe + m g- fc g^/2 G _ G+ gLo G G+ ^+_H g Wo/2 

= Q- k q w °/ 2 G„G + Q L °G„G + q w °' 2 vtJ:l 2k 
thus obtains the relations 

Vl k) g = 9Q- k vlr m % k . (8.31) 

(18.23P are contained therein as a special case where k — and m > 1. 
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When m = and k > 1, (I8.3ip become the relations 



V ik) g = gQ- k V^ ) . (8.32) 

According to the general rule [2TJ [221 EE], Vm^ and V^/c^ correspond to 
q-km/2qkM jjn an( ^ ^-km , '2 ^kM jm ^ an( j ^hese re l a tions imply the equations 

q kM = Q- k q- k2 L- 2k q~ kM , k = 1, 2, • - • , (8.33) 

for the Lax and Orlov-Schulman operators. Taking into account the twisted 
canonical relations (16.31) . one can easily see that these equations reduce to 
the single equation 

q M = (Qg)- 1 !- V 71 . (8.34) 

Thus the special solution of the Toda hierarchy determined by (I8.22p 
turns out to satisfy (18.271) and (I8.34j) . These equations may be thought of as 
generalized string equations in the Toda hierarchy. Lastly, setting q = e~ , 
rescaling M and M as (16.101) and letting H — > 0, these equations turn into 
the generalized string equations (18.11) for the dispersionless Toda hierarchy. 



9 Concluding remarks 

We have identified the Seiberg-Witten curves (I5.10p and (17.91) for the de- 
formed U(l) gauge theories. What about the Seiberg-Witten differential 
and prepotential? This issue is briefly discussed in our previous paper pQ 
(though not in a fully correct form, because of the incorrect formulation of 
the Riemann-Hilbert problem). We reconsider this issue in the present setup. 

After the idea proposed in the previous paper, we consider the critical 
values 

£ 4D(0) =£ 4D(0) [p (0) ]; £ 5D(0) = ^(0)^(0)] 

of the energy functionals as functions of (s, t). As the following calculations 
demonstrate, they may be thought of as the prepotentials. 

Let us show that the i-derivatives of these functions can be expressed as 

' 1 du^—p^\u), (9.1) 



dt k J uo k + l 



nc-5D(0)/ .x i* r u\ -Rku 

du—— p W'{u). (9.2) 



dt k J UQ -Rk 
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Since the treatment of the 4D and 5D cases are parallel, we illustrate the 
calculations for the 4D case only. Since p^ is a solution of (I4.12p . the 
identity 

PP f * dvg^"(\u - v\)p^'(v) + = («o < u < Ul ) 

holds. Integrated twice with respect to u, it turns into the identity 

j dvg^Qu - v\)pW(v) + ~ £ + d« + C 2 = 

(uo <u< ui), (9.3) 
where C\ and C*2 are functions of t only. Differentiating the explicit form 

J J uo 



of the critical value with respect to yields 



((.() 



^ 2 /:^(, M -^f>, +£ ^ 



+ / ^fn^'N- 



We can use (19. 3p to simplify the part inside the parenthesis as 

£4D(0) . ttl 9p (0)/ (u) , U1 „ fe+1 



<9f 



A: 



r j d P r{u),„ ^. p , fov . . 

/ cfa pi± L (C lU + C 2 )+ / dw— — pl 0) ' (u . 

y«o <^ Ju k + i 



On the other hand, since pl°^ satisfies (14. 3p . the identities 

fU\ 

dupf\u) = -\ \ duupW>(u) = 

'uo Juo 

hold. Differentiating them with respect to tk gives 

,< 0, '(«) „ r> a P i 01 ' 



"' A M_M = o, r duu M>> = o. 



10 The variational equations (|4.10p and (|4.25j) can be recovered from (|4.12l) and (|4.27|) 
by differentiating once. Therefore C\ can be identified with v/2. 
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Therefore the integral containing C\U + C2 vanishes, and we obtain (19. ip . 

By f)4.17p , we can rewrite the right hand side of ( 19. ip and (19. ip to contour 
integrals as 

dS 4D{0) (s t) 1 f z k+1 

at * [S,t) - 1 1 dz^—dWiz), (9.4) 



dtk 2ni J c k + 1 

dS r f i0) (s,t) 1 f , e~ Rkz 



r dz——dW(z), (9.5) 
dt k 2m J c -Rk K J K J 

where C is a simple closed curve encircling the interval [uq, Ui] anti-clockwise. 
Note that dW(z) is a single- valued differential, so that the contour integrals 
are meaningful in the usual sense. Unfortunately, W(z) is multi- valued on the 
cut plane, so that one cannot do integration by part to replace the differentials 

as 

~k+l C —Rkz 



dW(z) -z"W(z)dz, —^-dW(z) ->■ -e~ Kkz W(z)dz. 



k + 1 w w -Rk 

Apart from this problem, W(z)dz may be thought of as a candidate of the 
Seiberg-Witten differential. Actually, one can add the exact form dV(z)/2 
to replace it by 

— i-L = W(z)dz H -^- L , 



where S(z) is the 5-function ( 14.18)) . Thus we eventually reach the same 
conclusion as Marshakov and Nekrasov [T2"| IT3] that the total differential of 
the S function gives the Seiberg-Witten differential. The numerical factor 2 
is correlated with the denominator of (I6.34p and (I8.15P : if T^'s (or TVs) are 
used in place of t^s, this numerical factor disappears. 
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